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The curl (B) modes of cosmic microwave background (CMB) polarization anisotropics are a unique 
probe of the primordial background of inflationary gravitational waves (IGWs). Unfortunately, the 
B-mode polarization anisotropics generated by gravitational waves at recombination are confused 
with those generated by the mixing of gradient-mode (E-mode) and B-mode polarization anisotropies 
as CMB photons propagate through the Universe and are gravitationally lensed. We describe here 
a method for delensing CMB polarization anisotropies using observations of anisotropies in the 
cosmic 21-cm radiation emitted or absorbed by neutral hydrogen atoms at redshifts 10 to 200. 
While the detection of cosmic 21-cm anisotropies at high resolution is challenging, a combined 
study with a relatively low-resolution (but high-sensitivity) CMB polarization experiment could 
probe inflationary energy scales well below the Grand Unified Theory (GUT) scale of 10 16 GeV - 
constraining models with energy scales below 10 15 GeV (the detectable limit derived from CMB 
observations alone). The ultimate theoretical limit to the detectable inflationary energy scale via 
this method may be as low as 3 x 10 14 GeV. 
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Introduction — The curl (B) modes of cosmic mi- 
crowave background (CMB) polarization anisotropies are 
a unique probe of the primordial background of cos- 
mological gravitational waves 0. At these long wave- 
lengths, inflation 0] is the only known mechanism to 
causally generate such a background of gravitational 
waves. Since the amplitude of these inflationary grav- 
itational waves (IGWs) is proportional to V, the value 
of the inflaton potential V(^) during inflation, the am- 
plitude of gravitational-wave induced B-mode polariza- 
tion anisotropies directly constrains the energy scale of 
inflation V 1 / 4 (see, for example, Ref. 0). While the 
experimental sensitivity to B-mode polarization can be 
improved, the expected signal is contaminated by fore- 
ground effects |4| . The main confusion to the detection of 
B-mode polarization anisotropies generated by IGWs at 
recombination is the mixing of gradient-mode (E-mode) 
and B-mode anisotropies via gravitational lensing 5] . 

In this Letter we propose a new method for separat- 
ing lensing-induced B modes from the IGW signal using 
observations of anisotropies in the cosmic 21-cm radi- 
ation emitted or absorbed by neutral hydrogen atoms 
at redshifts 10 to 200. While the detection of cosmic 
21-cm anisotropies at high resolution is challenging, a 
combined study with CMB polarization data could probe 
inflationary energy scales well below the Grand Unified 
Theory (GUT) scale of 10 16 GeV — constraining infla- 
tionary models with energy scales below 10 15 GeV. The 
ultimate theoretical limit to the minimum detectable en- 
ergy scale of inflation via this method could reach as low 
as 3 x 10 14 GeV. 

Gravitational Lensing — Lensing induces a remapping of 
the polarization field at the last-scattering surface ±X(n) 
such that ±X(h) = ±X[n+V<f>(h)] is the observed polar- 



ization field, where ±X = Q±iU are linear combinations 
of the Stokes parameters Q and U and a(n) = V</>(n) is 
the lensing deflection angle. Here, 

0(n;z s ) = -2 / dr'^— $(n,r') (1) 
Jo r'r 

is the deflection potential, a line-of-sight projection of the 
gravitational potential $ to redshift z s . The total lensing 
potential </>(n) = 0(n; zcmb) is this quantity evaluated at 
z s — > ZCMB ~ 1100. 

Using the flat-sky approximation and the E-mode/B- 
mode decomposition [l|, the lensed B-mode polarization 
power spectrum, in the relevant limit Cf B << C EE , is 

qbb = c bb + J ^[l".lf sin 2 (2^C^, (2) 

where 1" = 1 — 1'. The second term in this expression 
is the lensing confusion in the B-mode map which must 
be separated from C BB — the IGW signal. Here, Cf 
is the angular power spectrum of the total deflection po- 
tential and is simply related to a weighted projection of 
the matter power spectrum 0; Cf^{z s ) is the incomplete 
power spectrum out to source redshift z s < zcmb- 

Unlike the B modes generated by tensor perturba- 
tions (the IGWs) , C EE is dominated by larger amplitude 
scalar perturbations. The expected few-percent conver- 
sion of E modes creates a large signal in the B-mode 
power spectrum [jj. For tensor-to-scalar ratios T /S be- 
low 2.6 x 10" 4 or, since V 1/4 = 3.0 x 10~ 3 (T "/ 'Sf'Snpi 
0, V 1/4 below 4.6 x 10 15 GeV the IGW signal is com- 
pletely confused by the lensing contaminant 0. To by- 
pass this limit one must separate the lensing induced 
B-modes from those due to IGWs. Clearly, the lens- 
ing confusion could be exactly removed if one knew the 
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three-dimensional distribution of mass out to the CMB 
last-scattering surface. However, as our goal is a mea- 
surement of C BB , knowledge of the projected quantity 
4>(n) is sufficient. One way to estimate 0(n; z s ) is by us- 
ing quadratic estimators or maximum likelihood methods 
to statistically infer the deflection-angle field given some 
lensed random field x(n) at source redshift z s . For in- 
stance, arcminute resolution CMB temperature and po- 
larization maps could be used to make such an estimate 
jil IflL fiHl . Another way to estimate 0(n; z s ) is by observ- 
ing the weak-lensing distortions of the shapes of objects 
of a known average shape at source redshift z s . We note 
here that observations of the weak lensing of galaxies, 
which have z s ~ 1-2, can not be used to delense CMB 
maps because a large fraction of the lensing contamina- 
tion (55% at I = 1000) comes from structure at z > 3. 
Higher source redshifts are required for effective delens- 
ing. We first review the potential observational signa- 
tures of the cosmic 21-cm radiation, and then discuss the 
methods for delensing the CMB outlined above. 
Cosmic 21-cm Radiation — Neutral atoms kinctically de- 
couple from the thermal bath of CMB pho tons at z ~ 200 
and cool adiabatically as T g oc (l+z) 2 [llj. Since the spin 
temperature of the hydrogen atoms remains collisionally 
coupled to T g these atoms resonantly absorb CMB pho- 
tons at A21 = 21.1 cm — the hyperfine transition of the 
ground state of hydrogen. The cosmic 21-cm radiation is 
thus first observable in absorption by low-frequency ra- 
dio telescopes which could detect brightness-tem pera ture 
fluctuations at wavelength A = A21XI + z) [lil ll3L H3 |. 
During reionization, the neutral gas distribution is likely 
to be complex due to the first luminous sources and 
cosmic 21-cm signatures shift to emission [l^. Yet, even 
before reionization, it is possible the 21-cm sky is bright- 
ened by emission from neutral hydrogen gas contained in 
minihalos with masses ~ 10 3 -10 7 M [13. 

Like the CMB, the statistics of the high-z absorption 
fluctuations are expected to be Gaussian and quadratic 
estimators of the lensing potential, described below, 
could be straightforwardly adapted to reconstruct the de- 
flection field. If the statistics of the 21-cm fluctuations 
during reionization can be understood the lensing of the 
21-cm emission from that era might also provide a useful 
probe of the lensing potential. The most promising (but 
futuristic) possibility involves using the shape statistics 
of high-rcdshift minihalos to infer the lensing potential. 
Quadratic Estimators — Quadratic estimators can 
be used to extract lensing information from the 
gravitationally-lensed field x(n) °f some intrinsic field 
x(n) at redshift z s . The quadratic form V • [x(n)Vx(n)] 
provides an estimate of the deflection angle at position n 
on the sky given the x anisotropy map. For the CMB, the 
quantity x could be the temperature anisotropies || , the 
polarization anisotropies, or some combination of both 
[lOJ . The brightness temperature fluctuations in the 21- 
cm transition of neutral hydrogen from redshifts 10 to 



200 could similarly be used |l8|. 

In Fourier space, the quadratic estimator for the de- 
flection potential is 
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where C xx is the unlensed power spectrum and 
T xx = C xx + N xx is to total power spectrum, includ- 
ing lensing corrections and a noise power spectrum N xx . 
The expectation value of the deflection-potential estima- 
tor (<j)(l;z s )} (the ensemble average over realizations of 
the random field x) is j us t 4>(\'i z s ). Here, 
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is the noise power spectrum associated with a quadratic 
reconstruction of Cf <! '{z s ) using the field x vk- 
Partial Delensing Bias — An estimate of 0(n) can be 
used to delense the CMB B-mode polarization map. In 
the limit z s — > zcmb (conventional CMB delensing) the 
extraction of C BB from the delensed map is limited by 
the noise introduced during delensing. The residual con- 
tamination of the B-modes is given by the second term 
of Eq. J5J) with the replacement Cf^ — > Qi . However if 
z s < zcmb this noise is not necessarily the factor limiting 
a measurement of the IGW signal. Using (f>(h; z s ) as a 
proxy for (j>(n) to delense the map leaves a residual lens- 
ing contamination not due to noise. Accounting for this 
partial delensing bias Bi(z s ) = Cf^ — Cf^ \z s ) (due to the 
difference in source redshift between the lensed field x(n) 
and the CMB) the residual contamination of the B-mode 
power spectrum is instead the second term of Eq. J2J with 
C; — > Bi{z s ) + Mi{z s ). This is true whether $(n; z s ) is 
estimated using quadratic estimators or by some other 
method. Here, Ni(z s ) is the residual noise power spec- 
trum of the deflection potential due to noise associated 
with the delensing process — for quadratic reconstruc- 
tion Mi{z s ) = Qi(z s ). If the deflection potential is re- 
constructed from a line, as is the case for cosmic 21-cm 
radiation, the source redshift is exactly know and Bi(z s ) 
can be reliably estimated. 

Quadratic Reconstruction — Unlike the CMB 
anisotropies, which lack power on angular scales 
below a few arcminutes due to Silk damping, the cosmic 
21-cm anisotropies extend to much higher values of I 
(limited by the Jeans wavelength of the gas) and peak 
in amplitude at higher values of I [HQ. Additionally, 
measurements of cosmic 21-cm anisotropies in different 
frequency bins provide several estimates of essentially 
the same deflection field. 

As shown in Fig. ^ we estimate a 21-cm experiment 
centered around z s ~ 30 with a 20 MHz coverage in fre- 
quency space capable of observing anisotropies out to 
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I ~ 5000 would have an Ni{z s ) higher than the planned 
CMBpol mission. In this case residual confusion arises 
from noise rather than bias. A quadratic reconstruction 
of the deflection field using this type of measurement, in 
conjunction with a CMB polarization experiment, could 
detect T/S > 2.5 x 10~ 5 or V 1/4 > 2.6 x 10 15 GeV. A 
next-generation 21-cm experiment capable of observing 
anisotropics out to I ~ 10 5 with the same central redshift 
and bandwidth would have a J\fi{z s ) an order of magni- 
tude below CMBpol and now be limited largely by the 
bias Bi(z s ). Paired with CMB polarization observations, 
this type of measurement could detect T/S > 1.0 x 10~ 6 
or V 1//4 > 1.1 x 10 15 GeV — comparable to very high 
sensitivity and resolution future CMB observations alone 
(lij ]. However, as the CMB data would not need to be 
used to reconstruct the deflection field a much lower res- 
olution CMB experiment would suffice. Furthermore, if 
lcnsing information need not be extracted from the CMB 
observations, the optimal observing strategy is to inte- 
grate over a few square degree patch of the sky as pro- 
posed in Ref. j2l|. Such a CMB experiment could thus 
be ground based. 

Other Methods — If minihalos bright in 21-cm emis- 
sion exist in the early Universe, just as galaxy shapes 
are sheared by weak gravitational lensing so will be the 
shapes of these minihalos. Ellipticity information ob- 
tained from such 21-cm minihalos could be used to recon- 
struct the projected potential out to high z [22^. Based on 
the dark-matter halo mass function, we expect roughly 
a surface density of 10 11 /sr of such minihalos at z ~ 30 
for a bandwidth of 1 MHz with masses between 10 5 and 
10 7 Mq. A typical halo of mass 10 6 M Q has a charac- 
teristic projected angular size of ~ 60 milliarcseconds. If 
resolved, then techniques currently applied to measure 
shear in background galaxies in the low-z Universe could 
be adapted for this application. Regardless of the exact 
method, a deflection potential reconstructed from high- 
redshift 21-cm observations could then be used to delense 
CMB B-mode maps. 

Bias-Limited Delensing — To understand to what extent 
bias-limited reconstructions, where the residual lensing 
contamination is dominated by Bi(z s ), would result in the 
removal of lensing confusion we have calculated the resid- 
ual B-mode power spectrum after correcting for the mod- 
ified lensing kernel when z s < zcmb- We have adapted 
the formalism of Ref. j2jj to estimate the smallest de- 
tectable background of IGWs and the resulting limits are 
summarized in Fig. 2. While knowing the projected mass 
distribution out to z s = 1 does not allow the confusion 
to be reduced significantly, if it is known to z s = 10 the 
confusion is reduced by an order of magnitude and the 
minimum detectable energy scale of inflation is reduced 
below the limit derived using quadratic CMB statistics 
[2oj | . A lensing-source redshift z s > 30 would be required 
to improve beyond the practical 1.1 x 10 15 GeV limit of 
the more sophisticated maximum likelihood method with 
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FIG. 1: Shown is angular power spectrum of the deflection 
potential as a function of source redshift z s . The curves la- 
beled 'I < 5000' and 'Z < 10 s ' are the estimated noise lev- 
els for quadratic reconstruction using 21-cm anisotropics in 
40 0.5 MHz bins centered around z s « 30. We assume a 
noise power spectrum with T sy3 = 3000iS' at 46 MHz, that 
Zmax/cov ~ 15, and a year of integration [T^. The curve 
labeled shapes shows the residual noise curve in a scenario 
where shear is directly measured using resolved minihalos in a 
1 MHz bandwidth about z 3 w 30. Also shown is the noise lev- 
els for a CMB reconstruction of deflections with the planned 
CMBpol mission assuming a 3 arcminute beam, a noise level 
of lfiK y/sec, and a year of integration. 

high resolution and sensitivity CMB polarization obser- 
vations |l9| . However we emphasize here that in the case 
where a bias-limited reconstruction of the deflection field 
exists for z s ~ 30, a very high-resolution CMB polariza- 
tion experiment is not necessary, and a similarly sensitive 
experiment with a lower resolution could do the same 
job. Thus, it is conceivable that the Inflationary Probe 
of NASA's Beyond Einstein Program can be designed in 
combination with a cosmic 21-cm radiation experiment. 
Such a cosmic 21-cm radiation experiment could be ex- 
tremely exciting in its own right 0, |2^|. For a bias- 
limited reconstruction out to a z s ~ 100, the limit on the 
tensor-to-scalar ratio is 7.0 x 10~ 8 or V 1/4 > 6.0 x 10 14 
GeV. For z s ~ 200, the maximum redshift where 21-cm 
fluctuations are expected to be nonzero, one could probe 
down to V 1 / 4 > 3 x 10 14 GeV. 

Discussion — While obtaining bias-limited measure- 
ments out to z s w 100 is a daunting task, with many 
experimental and theoretical obstacles to overcome, there 
is great interest in detecting the fluctuations in the cos- 
mic 21-cm radiation at z ~ 10 — 200 for their own sake. 
The observational study of 21-cm fluctuations, especially 
during and prior to the era of reionization, is now being 
pursued by a variety of low frequency radio interferom- 
eters such as the Primeval Structure Telescope (PAST 
[3), the Mileura Widefield Array (MWA), and the Low 
Frequency Array (LOFAR |2q). Planned interferometers 
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FIG. 2: Shown are power spectra of CMB B-mode polariza- 
tion. The curve labeled 'IGWs' is the IGW contribution to 
B-modes assuming a tensor-to-scalar ratio of 0.1 with (solid 
line; r = 0.17) and without (dashed line) reionization. The 
curve labeled 'lensing' is the total lensing confusion to B- 
modes. Thin lines show the residual B-mode lensing con- 
tamination for bias- limited delensing out to z s . Thick lines 
show previous estimates of the residual confusion from CMB 
experiments alone using quadratic estimators with an ideal 
noise- free experiment (dashed line), and likelihood methods 
using a high resolution/sensitivity experiment (dot-dashed). 
The noise curve of latter this experiment, with 2 arcminute 
beams and a pixel noise of 0.25 /^K-arcminute, is the curve la- 
beled 'High-Res'. Bias-limited lensing information to z s > 10 
improves upon the limit to the IGW amplitude based on 
quadratic statistics, and the likelihood level can be reached 
with z a ~ 30. If z s ~ 100 an additional order-of-magnitude in 
the IGW amplitude could be probed. The curve labeled 'Low- 
Res' is the noise curve for a lower resolution CMB polarization 
experiment with 30 arcminute beams sufficiently sensitive to 
detect IGWs when paired with a cosmic 21-cm lensing recon- 
struction. For very efficient delensing other foregrounds, such 
as patchy reionization |27|. might dominate confusion. 



such as the Square Kilometer Array (SKA) will improve 
both sensitivity and low- frequency coverage. 

While certain models of inflation, those related to 
Grand Unified Theories (GUTs), are expected to have 
an energy scale V 1/4 between 10 15 GeV and 10 16 GeV 
there are certainly other possibilities. For instance, some 
supersymmetric theories of inflation have energy scales of 
several times 10 14 GeV 2q| . New methods, such as the 
idea of using observations of the cosmic 21-cm radiation 
to delense the CMB B-mode polarization suggested here, 
are needed to push the minimum detectable energy scale 
of inflation below 10 15 GeV and discriminate between 
between physical theories at the highest energy scales. 
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